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Abstract 

The Q 2 evolution of polarised parton fragmentation functions is discussed using 
Altarelli-Parisi evolution equations. The first moments of both polarised quark and 
gluon fragmentation functions are shown to behave in a similar fashion at very high 
energies. This analysis is applicable to any hard processes involving production of 
polarised hadrons. The polarised A hyperon production in e + e~ annihilation where 
this can be realised is considered. We present a complete a s (Q 2 ) corrections to the 
asymmetries discussed in the paper of Burkardt and Jaffe which demonstrates the 
extraction of various polarised fragmentation functions. To this order, these corrections 
are found to be scheme dependent similar to that of structure functions. 



1 Introduction 



The theory of strong interaction physics is well understood within perturbative Quantum 
Chromodynamics(PQCD) [|IJ. The PQCD describes the dynamics of quarks and gluons at 
very high energies. The laboratory experiments involving hadrons either in the initial state 
or final state can be understood by treating these quarks and gluons as asymptotic states and 
by safely using the perturbative techniques thanks to asymptotic freedom. In other words, 
the hadronic cross sections can be well expressed in terms of partonic cross sections. The 
Parton Model (PM) [Q] has been a successful model to tie up these perturbatively calculable 
parton level cross sections and the experimentally measured cross sections. This involves the 
introduction of certain probability distributions which can not be calculable in PQCD but 
are just inputs of the model. Though they are not calculable, they are found to be universal 
in the sense that they are process independent and the evolution in terms of the scale is well 
known. In other words, if these probability distributions are measured in an experiment at a 
scale, say Q 2 which is above the Aqcd where the QCD perturbation theory is reliable, then 
these distributions can be used as inputs to predict the rates or cross sections of different 
experiments which involve these distributions. In this sense they are something to do with 
the hadrons participating in the process and are process independent. Also, using the well 
known Altarelli-Parisi evolution equation satisfied by these distributions, one can find out 
how these distributions change as the scale changes and thereby use them appropriately for 
the experiments done at various energies. 

There are two types of such probability distributions one comes across in hadron physics. 
The well understood among them is the parton probability distribution. It is generally de- 
noted by f a (h)/H(s)( x ,Q 2 ) an d the meaning of it is the probability of finding a parton of 
type a with polarisation h inside the hadron H of polarisation s with momentum fraction 
x of the hadron and the scale Q 2 . These distributions are usually measured in Deep Inelas- 
tic Scattering(DIS) experiments. The universality of these distributions and the predicted 
evolution in terms of the scale Q 2 are proved experimentally . The other distribution is 
nothing but the fragmentation function. The fragmentation functions are mirror image of 
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the parton probability distributions. These functions are generally denoted by D a ^J(x, Q 2 ). 
They measure the probability that a parton of type a with polarisation h fragments into a 
hadron of type H with polarisation s and carries away x fraction of the parton's momentum 
at a scale Q 2 ||. These functions are usually measured in e~ e + — > H X experiments. 
The universality and the Q 2 evolution of unpolarised functions (which does not contain any 
information about the spin of the hadrons produced) are well understood. 

In recent years, there has been several interesting works to understand these fragmenta- 
tion functions. These include the measurement of the fragmentation functions for charged 
and neutral pions and kaons The QCD inspired Parton Model analysis sublamented 
with Altarelli-Parisi evolution equations for these functions is shown to explain the pion and 
kaon production rates at leading and next-to- leading order ||. The another field of current 
interest is the inclusive A production at e~ e + scattering experiments || . This is important 
as A s are produced copiously and also it is easy to detect them. On the theoretical side, 
the study of polarisation effects in the production cross sections is more interesting. As 
we know, at low energies the unpolarised e~ e + scattering will not produce longitudinally 
polarised A s (we do not attempt to review the transversely polarised A s produced in unpo- 
larised scattering experiments) but when the center of mass energy crosses the Z threshold, 
then due to parity violating effects, the produced A s can naturally be polarised. Recall that 
these are produced due to the fragmentation of polarised partons produced in the scattering. 
Unfortunately, it is not clear how these partons fragment into these hadrons. There exists 
only model calculations which can tell us how much of the parton spin is transferred to A. In 
the naive quark model, the only contribution to polarised A is due to strange quark and the 
contribution coming from other partons is identically zero. As we know this naive picture 
breaks down due to complicated QCD effects of both perturbative and non-perturbative 
origins. Along this line of thought there has been an interesting work by Burkardt and Jaffe 
[0] who have charted out an experimental programme to measure polarised u and d quark 
fragmentation functions in addition to s quark fragmentation function. The non-zero value 
for u and d quark fragmentation functions will invalidate the quark model picture. More 
recently the references || , || discuss the importance of gluons in the polarised A production. 
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The analysis in the reference f8[ is purely based on the AP evolution equations satisfied by 
the polarised fragmentation functions of quarks and gluons. It has been shown that the 
gluons play significant role when the scale at which the experiment is done is very high. In 
the massive gluon scheme, it has been demonstrated || that the gluons contribute to po- 
larised lambda production. Since the fragmentation functions are defined in this scheme, the 
gluonic contribution to polarised A is scheme dependent. The situation is very similar to the 
gluonic contribution to polarised structure function measured in DIS. This paper extends 
the previous analysis to include the Z boson exchange in order to completely extract 
various parton fragmentation functions. 

The paper is organised as follows. In the second chapter we discuss the importance of 
polarised gluons in the A production using AP evolution equations. We then systematically 
compute the QCD corrections to the asymmetries defined below which are useful to extract 
various partonic contributions to the A production. The results are presented in the chapter 
3. We finally conclude in the chapter 4. The appendix contains the relevant details of the 
results presented in the chapter 3. 

2 The Q 2 Evolution of the Polarised Fragmentation 
Functions 

In the parton model, the production cross section for A with polarisation s in e~ e + scattering 
is related to that of a parton and the probability that the parton fragments into A. In other 
words, 



The left hand side is the hadronic cross section for the production of polarised A within a 
solid angle Q with respect to the beam direction in e~ e + scattering. Here s, s' and A are the 
helicities of A, electron and parton respectively. The Bjorken like variables x = Ip^.qjQ 2 
and y = 2p a .q/Q 2 are scaling variables in the hadron and parton levels respectively with 
q 2 = Q 2 , pA,p a and q being the momenta of produced hadron, parton and the intermediate 




(1) 
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vector boson(7 or Z) respectively. The production cross section for a parton with helicity A 
(right hand side of the equation) is completely calculable in PQCD. On the other hand, the 
probability that a parton of helicity A fragments into A of helicity s, D^k}(x/y,Q 2 ), is not 
calculable due to non-perturbative mechanism involved in the fragmentation region. Hence, 
these fragmentation functions are just inputs of the model which can be extracted from the 
experiment. Though they are not calculable within the realm of PQCD, they are universal 
in the sense that they are process independent. More interestingly their evolution in terms 
of Q 2 and x are completely governed by the well known Altarelli-Parisi evolution equations. 

In the polarised experiments, the hadronic cross section is directly proportional to po- 
larised parton fragmentation functions defined as 



AZ^x, Q 2 ) = D$ (x, Q 2 ) - D*l (x, Q 2 ) 



(2) 



Let us first analyse the Q 2 evolution of these fragmentation functions using Altarelli- 
Parisi evolution equations. The evolution equations are well known and are given by 



dt gi 2ir Jx y 

at y 2n Jx y 



AD$(y,t)AP qq (x/y) + AD$(y,t)AP gq (x/y) 



2/ 



£Al£ (y, t)AP qg (x/y)+AD$(y, t)AP gg {x/y) 

3=1 



(3) 



(4) 



where t = log(Q 2 /A 2 ) and a s (t) is the strong coupling constant. Here, a s (t)AP ab (y)dt is the 
probability density of finding a parton of type a at the scale t + dt with momentum fraction 
y inside the parton of type b at a scale t. The splitting functions are given by 

AP qq (z) = c 2 (R)(^±^- + p(l-z)\ 



C 2 (R) 
1 



z 2 -(l- z) 2 ) 



AP gg {z) = C 2 (G) \{l + z 



+ 



>\z (1-4 

11 _ I H^L\ 

~6~3C 2 {G)J 



6(1 -z) 



Here, C 2 (R) = (N 2 - 1)/2N, C 2 (G) = N and T(R) = f /2 with A^ = 3 for 577(3) and / being 
the number of flavours |IJ]. In the above equations, usual + prescription has been used to 
regulate z — > 1 singularity. Notice that the above equations are similar to the AP equations 
satisfied by the polarised parton distribution functions but for the interchange in Ap qg and 
Ap gq . The reason for this is as follows: The emission of a quark(gluon) from a quark (gluon) 
only affects the probability of quark (gluon) fragmenting into hadron. Hence the splitting 
functions AP qq and AP 59 are unaffected. On the other hand, the emission of a quark from a 
gluon changes the probability of the gluon fragmenting into hadron. Similarly, the emission 
of a gluon from a quark affects the probability of quark fragmenting into hadron. That is 
why, the splitting functions AP qg and AP OT are interchanged. As these equations can easily 
be solved in the Mellin space, we define 

AD A (n,t) = [ x n - l AD%(x,t)dx 
Jo 

AP ab (n) = f x^APatWdx (5) 
Jo 

The complete solution for the nth moment is not illuminating. Recall that the first moments 
of the polarised parton distributions are interesting as they are related to the spin content 
of the polarised hadron and the measurement of them will tell us how the hadron spin is 
shared among the partons. In the same spirit, we here look at the first moment of polarised 
parton fragmentation functions. The measurement of it for various hadrons will tell us how 
the parton helicity is distributed among the produced hadrons. But it is experimentally a 
hard task. From eqn.([|), we find that the first moment of the polarised gluon fragmentation 
function to order a s (t) satisfies a simple first order differential equation, that is 

^AD^t)=a s (t)(3 AD^t) (6) 

where fi = (11C 2 (G) — 4T(P))/127r. The solution to the above equation can be found very 
easily using renormalisation group(RG) equation for the QCD coupling constant, 

ja s {t) = -(3 a 2 s (t) (7) 

From eqns.(|| 0), we obtain an interesting behaviour of first moment of gluon fragmentation 
function: the product of the first moment of polarised gluon fragmentation function times 
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the strong coupling constant is scale independent to order a 2 (t), 

j t (a s (t)ADj(t)) = 0(a s (tf) (8) 
In other words, to order a 2 (t), ADg increases as the scale t increases, i.e 

ADi(t) = K]og(9pJ (9) 

where K is some constant. It is worth recalling that the counter part of such a relation for 
polarised gluon distribution function exists and has opened up a better understanding of the 



spin structure of the proton [10|. That is, 

j t (a s (t)Ag(t)) = 0(a s (t) 2 ) (10) 

where Ag(t) is the first moment of polarised gluon distribution function. From the above 
equation it is clear that polarised gluonic contribution to proton spin could be important 
at very high energies. But this equation does not say anything about the absolute value of 
gluonic contribution. 

Now let us consider the first moment of the polarised quark fragmentation function into 
polarised hadron. From eqn.(D), it turns out 

J t ^ D ^) = l<t)ADi(t) (11) 

From eqns. @ and (p]), we find that AD^.(t) grows as t. This is to be compared with the 
first moment of the polarised quarks in polarised hadron which is scale independent. 

Hence the relations(eqns.(^i|TTl)) suggest that at very high energies the polarised gluon 
fragmentation into polarised hadrons(say A) is as significant as polarised quark fragmen- 
tation. This point is very important when one is interested in the analysis of integrated 
asymmetries coming from various partons. Also the first moment of the fragmentation func- 
tions will tell us given a parton of definite polarisation, how much of its polarisation is 
transferred into the produced hadron. This is to be compared with the interpretation of the 
first moment of polarised parton distributions that it measures how much of hadron's(say 
proton's) helicity is shared by the parton. If we sum over all these contributions, it will turn 
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out to be the helicity of the hadron. Similarly in the fragmentation sector if we sum over 
all the hadron fragmentation functions coming from a specific polarised parton with their 
appropriate polarisations then it will coincide with the polarisation of the parton. 



3 Asymmetries 

The analysis in the last chapter shows that the gluonic contribution to the polarised A is 
as important as the quark contribution. Hence the QCD corrections to the processes which 
involve the extraction of these fragmentation functions are important. Recently, Burkardt 
and Jaffe J7| have discussed the method of extracting various polarised quark fragmentations 
to polarised A by measuring some specific asymmetries in both unpolarised and polarised 
e~ e + scattering experiments. As a preliminary effort, in the reference ||, |9|] various QCD 
corrections are computed to the cross section asymmetry when there is no Z exchange 
(purely electromagnetic (em)). The factorisation of IR singularities has also been verified to 
order a s (Q 2 ) |§. In the following, we compute the QCD corrections to those asymmetries 
(discussed by Burkardt and Jaffe) which involve both 7 and Z vector bosons as intermediate 
particles. 



3.1 The Unpolarised e e + Scattering 

First we compute the following asymmetry 0. 

- d& a(l) dd a(V 



da A ^ da A ^ /-i dy 



dxdVt dxdVt „ J x y 



AD A (x/y,Q 2 ) (12) 



dydVt dydVt 

In the above asymmetry the spins of the initial leptons are averaged. This asymmetry is 
identically zero when there is only 7 exchange. At very high energies, the Z exchange is 
also possible. This will make the asymmetry non-zero thanks to parity violation. In the 
above equation Q is the solid angle within which the produced parton fragments into A. The 
arrows in the parenthesis of A and the partons denote their helicities with respect to the 
beam direction. The scale Q 2 is the invariant mass of the photon or Z produced(i.e.,g 2 = Q 2 ) 
at the e + e~ vertex. The sum is taken over all the partons such as quarks, anti-quarks and 



gluons fragmenting into A. Recall that the kinematic scaling variables x and y are defined 
as x = 2p\.q/Q 2 and y = 2p a .q/Q 2 respectively. The parton level asymmetry given in the 
eqn.(P^D is factorisable as 

d&a (i) ^a(T) 
dydVL dydVL AQ jJ^ lZ 



^ E C^(Q 2 )V {I) (Q 2 ) E TfHf>{y,Q*) 



(13) 



3=1,3,4 



where V^(Q 2 ) are propagators given by 



Mi 



(14) 



(Q 2 -M 2 Y Q\Q 2 - 1V1Z) 

The tensors T- IV are constructed by looking at the symmetry properties of the amplitudes 
for direct Z and interference contributions. These two amplitudes are individually 
gauge invariant in the mass less limit as can be checked explicitly. Hence these tensors 
are parametrised as 



—e^qxPaa 



Pa-q 



2 3 



2 4 



1 



,Pa-Q 



(15) 



(16) 



Pa ■ q V - Q J V Q J 

The tensors proportional to q^ in T 3 and T 4 are immaterial as the leptonic tensors are 

individually gauge invariant. The leptonic tensors C^s can be easily workout and are found 



to be 



Hz) 



Sin 2 6 w Cos 2 8\y 
4ira 

SinOwCosOw 



(v 2 + a 2 X v ~ 2iv e a e F u 



(17) 



where 

l»» = q»q 2 ' + q»q 2 l -g^q 1 .q 2 

with qi and q 2 being the momenta of incoming positron and electron respectively. Here a is 
the fine structure constant, v e and a e are the vector and axial vector couplings in the e + e~ Z 
vertex and 9w is the Weignberg angle [IT]. It is clear from the above tensors that the terms 
proportional to q^ or q u in the tensors T 3 and T 4 when contracted with the leptonic tensors 
give zero contribution. 
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Substituting the tensors (eqn.fllBD) and the leptonic tensors (eqn. ([T7|) ) in the eqn.([T 
we obtain 

da a ^ da a{ -V 1 na 1 



CosOHl + H% + |(1 - Cos 2 9)Ht 



(19) 



«i = ^ m o^n —n v e a e H[ Z)a + 2Qa e Hf Z)a (21) 



dydQ dydfl 2Q 2 Sin9 w Cos9 w (Q 2 - Mf 
where the general form of TCf for i = 3,4 is given by 

U a = — - (v 2 + a 2 )H {z)a + 20v H {lZ)a (2G\ 

and that of Tif is given by 

ci 1_ 

Q 2 - M\ Sin9 w Cos6 w ' 

where, v q and a q are vector and axial vector couplings in qqZ vertex jOj] . The superscript in 
H denotes the origins of these contributions such viz, direct Z or the interference between 
7 and Z contributions. The contribution to H^ a can come from processes to lowest order 
(see fig.l)(a^°' ) ) in a s (no gluon emission) as well as from processes involving single gluon 
emission(see fig. 2) and virtual contributions(see fig. 3) to first order in a s . The evaluation of 
asymmetries to lowest order is very simple as it involves only two body phase space. When 
the gluon accompanies the quark and anti-quark pair, then the evaluation is cumbersome 
due to three body phase space integral. We have given below a simple looking formula to 
compute the three body phase space after performing most of the integrals using the delta 
functions. The H^' a s are related to the matrix elements as follows: 

#f)° = [ dxiVt u — \Mj\ 2 (i - t) (22) 

32(2?r) 3 i 1 4tt IIabA+ u ( j 

where the projectors are given by 



1 2p a .q 

y V Q y 



?T = r(^ + 12 ^l ( 23 ) 



The terms H\ I)a can be computed from the fig. 2 and are found to be of the form 



1 ~ 6 64TiSin 2 6 w Cos 2 6 w Ll ' ~ 6 167iSind w Cosd w ^ ^> 



where 

C\ Z)a = Wq Cf C^ Z)a = Va a q Cf 

with i = 3,4 and r\ q = 1 for quarks and —1 for anti-quarks. The factor 3 appearing in the 



eqn.(p4|) is due to the number of colours. The functions Cf are computed in the appendix. 
As we have already mentioned the functions Cf to lowest order get contributions from the 
Z or the jZ decaying into a polarised quark and antiquark pair. This asymmetry is just 
proportional to 5(1 — y) but for charge and other group factors. To order a s , Cf s get 
contributions from two types of processes: 1. polarised quark or antiquark production with 
a real gluon emission and virtual gluon corrections to the quark and anti-quark(self energy 
and vertex corrections) 2. polarised gluon emission from unpolarised quark and anti-quark 
pair. These processes suffer Infrared(IR) divergences when masses of the quarks and gluons 
are taken to be zero. In order to regulate these divergences we give gluons a small mass m g 
and finally take the limit m g going to zero. We compute these cross sections in the limit 
p a .q and Q 2 tending to infinity with their ratio fixed. This is analogous to DIS limit but the 
scaling variable here is inverse of the Bjorken scaling variable. 



4a/ 4 N 

C a 3 = W(1-^) + 3^(Q-^p 



C\ = -26 a , q/q 6(l-y) + ~^(--CZ + ^C a P 

3 2tt \ y y s 



4 (T 

C? = S a , q/q 6(l-y) + -^C- (26) 

Notice that the delta function is absent for gluons. The functions Cf for quarks are given 
below, 

d = d + 2-y 



C q P = \ (27) 



2 
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Notice that the above expressions are well defined in soft limit. As is shown in the appendix, 
the soft singularities coming from real gluon emission diagrams are exactly cancelled by 
that coming from the virtual diagrams such as self energy and vertex corrections. The term 
log(mjj) purely comes from the collinear divergence which can not be avoided in the mass 
less limit. This ill defined term will finally be absorbed into the fragmentation functions at 
the level of A production cross section. Hence the fragmentation functions are defined in 
what is usually called the massive gluon scheme. Similarly for gluons we have 

G{ = 2 (2-y)log(^f]-A(2-y) 



C 9 G = 



C 9 P = (28) 

(I)a 



The asymmetry given in the eqn.(|T2"D is expressed in terms of the C\ s as 



dxdVt dxdQ 



where 



E 



~-q,q,9 



x 



x 



+2xa(vi + ai)v q a qVq I C% + -(1 - Cos 2 6)C2 
-2Cos6 (xie q v q a e C± - X2deVe{v 2 q +a 2 q )Cl 



AD„(x) 



(29) 



Xi 



Q 2 



l 



Q 2 - Ml 16Sin 2 9 w Cos 2 9 



X2 



Q 4 



l 



(30) 



, z , 'iii , u, u W (Q 2 - M 2 Z ) 2 256SmW w Cos 4 6 w 

The Xi an d X2 for Z propagator with finite width will be modified and are found in ref. |11 
and the convolution of two functions f(x) ® g(x) means 



f(x)®g(x) = / —f{y)g{x/y) 
J x y 



(31) 
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3.2 The Polarised e e + Scattering 



Now let us find out the following asymmetry 
produced parton are polarised 



where one of the initial leptons and the 



dxdQ 



dxdQ 



E 



l dy 
y 



AD A (x/y,Q 2 



(32) 



dydVt dydfl 

Following the similar procedure, we can decompose the parton level asymmetry in terms of 
leptonic(£ /i!/ ) and gauge invariant hadronic(T All/ ) tensors as 



da a ^(i) da a ^(i) 



1 



(33) 



dydQ dydQ 

where the new V^(Q 2 ) = l/Q 4 . In the above equation, I = 7 corresponds to the diagrams 
where only photons are intermediate vector bosons. This is the extra contribution which 
is absent in the previous asymmetry. This is possible because of the polarisation of initial 
lepton. This diagram is gauge invariant as it is. The diagrams for I = Z, ^Z are discussed 
already. Notice that the diagrams are split in such a way that for each / the diagrams are 
individually gauge invariant, hence the simple tensor decomposition. The leptonic tensors 
^fn(i) f° r ^ ne polarised electron are found to be 



'(7) 

'(Z) 



2Sin 2 e w Cos 2 6 w 



[v e + a e 



2na 



(34) 



'{nZ)\+> Sin6 w Cos6 w 
where and are already given. Substituting the above leptonic tensors in eqn. (|33|) we 



find, 



da'Mjl) da^\j) _ 1 txol (v e + a e ) 

dydQ dydQ ~ 2Q 2 Sin6 w Cos6 w Q 2 - M 2 Z 



[HtiD+CosenKi) 



y -(i 

4 V 



Cos 2 9)Hta) 



(35) 



The general form of Tl^(i) is found to be 



Q 2 - M\ 2Sin6 w Cos6 w 



[v e + a e )H\ z)a + 2QH, 



r ( 7 Z)a 



(36) 
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The functions H^ a for i = 1, 3, 4 are given in the eqns.(p4|). The new 7i"(|) is given by 

= Air, (37) 



•57T 



Substituting eqn.(^) in eqn. (|32|) and using the eqns. (p6| , |37| , p4] , p5D we find 



dxdfl dxdVt 2Q< 



a r fx \ 

3 ^E [-Xi(ve+a e )a g e q r ]q [2C a 3 +-(l-Cos 2 e)C a 4 j 



+ X2 (v e + a e ) 2 v q a g r, q (2C% + |(1 - Cos 2 d)Ct 
-2 xi (v e + a e )v q e q CosO C" 
+X2 (v e + a e ) 2 (v 2 q +a 2 q ) CosO C\ 

+ e\ CosO Cf\ ® AD A {x) (3J 



From the eqns.(p9|, |38|) we find that the asymmetries reduce to those given Burkardt and 
Jaffe(BJ)'s paper when we put the strong coupling constant a s {Q 2 ) equal to zero confirm- 
ing the correctness of our results presented in this paper. Also, we can reproduce the results 
for the asymmetries in the absence of Z exchange diagrams by putting both \\ an d X2 equal 
to zero. Hence our results presented here are consistent with BJ asymmetries. The extrac- 
tion of both quark and gluons now becomes complicated because of the presense of gluonic 
fragmentation functions. Notice that the expressions (eqns. (p9|), (|38D) look very similar to 
that of BJ asymmetries with the replacement of quark and antiquark fragmentation func- 
tions by the parton fragmentation functions convoluted with complicated looking functions 
Cf. This simple fact that the structures are identical might help us in the extraction of the 
combination of fragmentation functions with Cf instead of just the fragmentation functions. 
From this combination and an independent measurement on polarised gluonic contribution 
one should be able to disentangle both the quark and gluon fragmentation functions. Once 
these functions are measured in the laboratory, these have nice physical interpretation that 
their first moments will tell us the polarisation contribution coming from various partons to 
a specific hadron in the sense described earlier. As is seen from the analysis of Q 2 evolutions 
of these fragmentation functions using AP evolution equations, the first moment of quark 
and gluon fragmentation functions will play a crucial role at high energies in the production 
of polarised A hyperons at e~ e + annihilation processes. 
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4 Conclusion 



In this paper we have extensively studied the polarised A production in e + e~ annihilation 
process. This involves the measurement of what are called polarised fragmentation functions 
of quarks and gluons into polarised A. The extraction of the fragmentation functions of 
different flavours is very complicated due to the charge factors multipling them. Burkardt and 
Jaffe were successful in disentangling these distributions by looking at various asymmetries at 
different kinematical regions. In this paper, using the Altarelli-Parisi evolution equation, we 
have shown that the polarised gluon fragmentation is also important when one is interested 
in finding spin or helicity coming from various partons to the produced polarised A. In 
fact we find that the first moment of the polarised gluon fragmentation is logarithmically 
raising as Q 2 increases. So at high energies gluons play crucial role invalidating the naive 
expectations based on simple coupling constant arguments. This is similar to the behaviour 
one encounters in the study of spin content of the proton or any hadron. In fact the behaviour 
in the fragmentation sector is much more stronger than that in the structure function sector 
due to the vanishing of first moment of AP qg (x). Since the gluons can contribute only from 
order a s (Q 2 ), we have considered the full QCD corrections to the asymmetries discussed 
earlier. This of course, further complicates the extraction procedure. The phenomenological 
implications of these QCD corrections and the extraction of gluonic fragmentation function 
are under investigation. It is worth to recollect the status of gluonic contribution to polarised 
structure function measured in DIS experiments. The status here is very similar to that 
in DIS due to the scheme dependence of these gluonic contributions. Hence, a scheme 
independent measurement of polarised gluonic fragmentation function would substantiate 
our analysis. The extraction of these distribution is very important primarily due to two 
reasons. One of them is that this might tell more about quark and gluon dynamics in the 
fragmentation sector. Particularly, the measurement of non strange quark contributions to 
polarised A production will be a test of existing models for polarised A fragmentation. Also, 
the importance of polarised gluons in the A production can be experimentally tested. 

I would like to thank prof. M.V.N. Murthy for his constant encouragement. It is a 
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A Matrix Elements and Integrals 

In this appendix we outline the evaluation H^ a in the massive gluon scheme. We have split 
the amplitude into three parts. That is, 1. 7 — > qqg, 2. Z — > qqg and 3. the interference of 
7 and Z initiated processes. All three are independently gauge invariant. 

Let us first consider 7 — > qqg process. If the quark is polarised then there will be virtual 
contributions due to gluons in addition to gluon emission diagrams. These diagrams are both 
soft and collinear divergent. In addition, the virtual diagrams are ultraviolet divergent (UV). 
The soft and collinear divergences are regulated by giving a small mass to the gluon. The 
UV divergence is regulated using Pauli-Villars regularisation. 

To produce polarised quark, unpolarised anti-quark and gluon the initial virtual pho- 
ton should be polarised, otherwise the amplitude is identically zero. We project out this 
contribution as 

ie» vXaP ^\M% v {l - t) = 2Waa s eyM\ (39) 
2p 2 .q 

where the momenta of the particles are given in the fig. 2, the Mandelstam variables are 
s — (pi + P3) 2 , t = (p 2 + P3) 2 and u — (pi+ p 2 ) 2 satisfying s + t + u = Q 2 + m 2 . Here, 



M\ = (st-m 2 g Q 2 )( 



s + 2t- 2m 2 g - Q 2 



+ ^7 



t 2 (s-Q 2 ) s\ 
+2(s + t-m 2 - Q 2 y - ' " 2 " 9 \ — ( 10) 



Q 2 s + st- m 2 g Q 2 - Q 4 
'•" " > st{Q 2 -s) 



The processes involving only Z particles are given in the fig. 2. Due to the axial coupling 
in addition to the vector coupling, polarised particles can be produced from the virtual Z 
decay even if one sums the polarisations of Z unlike the case of 7. Hence we project out the 
following contributions from the matrix element when the quark is polarised: 
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^D*" " T) = cXZeJ * + 4)^ (4i) 

a<& = (rt - m ;o») (^ + i) + 2( n ,; + <y) m ' + < ^— ' 

Here .M? is already given in the eqn.(pO|). 

Next we compute the interference term coming from Z and 7 decaying into polarised 
quark, anti-quark and gluon 



9 |A VUU \) Co80 w SinB w a * Ma 



2p 2 .q 1 * Cos0 w bm9 w 

Noting that s = Q 2 (l — x 2 ) and t = Q 2 (l — x\) where Xi = 2pi.q/Q 2 the relevant integrals 
are presented below 



/ 



dx x M\ = ( ] + X2 \ ( log(x 2 (l - x 2 ) + log ( ^ 

1 — X 2 / y y 7TI 

— (1 - x 2 ) - - — h -5(1 - x 2 ) 

2 V ; 2 1 -x 2 4 V y 



9- 



y dxiMc = Jdx\Mi + 2 — x 2 



dxiMp = X A (44) 

The virtual contributions can be computed from the diagrams given in the fig. 3. Here we 
state the expressions for the amplitudes involving self energy correction and vertex correction 
(for general vertex of the type (0 + 675)7^). Both diagrams(self energy and vertex correction) 
in fig. 3 are UV divergent separately. As we have already mentioned, we regulate these UV 
divergences in Pauli-Villars regularisation scheme. They also suffer from soft divergences 
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when we set all the masses to be zero. Here also, we give a small mass to the gluon to 
regulate them. The amplitude for the self energy insertion after properly taking care of wave 
function renormalisation turns out to be (self energies for both quark and antiquark) 



M self 

A* 



a. 



log 



m r 



v s ,(pi)(a + 675)7^(^2 



(45) 



where L is the Pauli-Villars regulator to regulate the UV divergence, s and s' are the 
polarisations of quark and antiquark respectively. 

Similarly the amplitude for the vertex correction turns out to be 



M, 



vertex 



+ - 



71 

7 



log 

3 



2 f m T 



Q 2 



I m n 

3l0g W 



log 



mi 



(46) 



Though the above amplitudes suffer from UV divergence, when we sum these two amplitudes, 
we find that the UV regulator cancels between self energy and vertex corrections, thanks to 
Ward identity. 

Now let us find out the projected matrix elements similar to eqns. (|40|) , (|42| ) when the 
gluon is polarised and quark and antiquarks are unpolarised. We here list Aif analogous to 
M.1 (with P2 replaced by p$) 



7] st- m 2 g Q 2 \ n (s + t)(s + t- m] - Q 2 ) ( Am 2 g Q 2 - 2m 2 g s - 2Q 2 s + s 2 - t 2 



st 



+ 



t 2 



4m 2 Q 2 - 2m 2 t - 2Q 2 t + t 2 -s< 



M 9 G = 


7] 1 


2m 2 g (s + t 


2s + t 


M 9 P = 




(mj + Q 2 - 


2Q 2 


s + t 



2m 



—t + (2Q 2 -s-t)(st-m 2 g Q 2 ) y ^ 

2 t- s 2m 2 Q 2 - m 2 s + m 2 t - 2st 
~t t 2 



1 

T 2 



2m 2 Q 2 + m 2 s — m 2 t — 2st 



(47) 



Where rj = 1 for 7 or Z diagrams and 77 = 2 for 7Z interference terms. We find that the 
results given in eqn. (^) can be got by integrating the above matrix elements with respect 
to %2- The integrations of the matrix elements M. 9 G and Ai 9 P are identically zero in the limit 



m g taken to be zero. 
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Figure Captions: 

1. Graph contributing to e + (qi)e~(q 2 ) — > q{pi)q{p2)- 

2. Graphs contributing to e + (gi)e"(g 2 ) — > Q.{pi)q.{P2)g{p?)- 

3. Self energy and vertex corrections to e + (g!)e _ (g 2 ) — * ?(pi)<z(p2) 



19 



